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ABSTRACT 

Context. The full three-year Wilkinson Microwave Anisotropy Probe results (hereafter WMAP3) reinforce the absence of large- 
angle correlations at scales greater than 60°. The Poincare dodecahedral space model, which may naturally explain such features, thus 
remains a plausible cosmological model, despite recent controversy about whether matched circle searches would or would not push 
the topology beyond the horizon. 

Aims. We have used new eigenmode calculations of the dodecahedral space to predict the cosmic microwave background (CMB) 
temperature anisotropics in such a model, with an improved angular resolution. 

Methods. We have simulated CMB maps and exhibited the expected self intersection of the last scattering surface along six pairs 
of circles. For a set of plausible cosmological parameters, we have derived the angular power spectrum of the CMB up to large 
wavenumbers. 

Results. Comparison of the angular power spectrum with the WMAP3 observations leads to an optimal fit with the Poincare dodeca- 
hedral space model, for a value Cli^t = 1.018 of the total energy density parameter. 
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' 1. Introduction 

. After lying dormant during the mid twentieth century, interest in 
• cosmic topology re-awakened with the COBE satellite observa- 
tions (Hinshaw et al. 1996 ) followed by the first-year Wilkinson 
Microwave Anisotropy Probe observations (Bennett et al. 120031 
hereafter WMAPl), which gave unexpected discrepancies with 
a ACDM model with infinite flat spatial sections (the so-called 

■ "concordance model"), including weak multipoles at ^ = 2 and 

■ 3, and violations of statistical isotropy in the same multipoles 
(Schwarz et al. 120041 ). Among various possible explanations, it 
has been suggested that a spatially finite multi-connected uni- 

' verse may explain the weak large-angle correlations (Hinshaw et 
al. 119961 Spergel et al. 2003). More precisely, when the spatial 
sections are taken to be the Poincare dodecahedral space (PDS), 
the low-^ cosmic microwave background (CMB) multipoles fit 
well with the observational data (Luminet et al. 2003), although 
it has been stated that the PDS model does not account for the vi- 
olations of statistical anisotropy (Aurich et al. 2007 ). More gen- 
erally, it was proved that the long-wavelength modes are rela- 
tively lowered in so-called "well-proportioned spaces" (Weeks 
et al. 12004 1). Along with the PDS model, the well-proportioned 
spaces include roughly cubical flat tori and the binary tetrahedral 
and octahedral spaces. The power spectra for the octahedral and 
tetrahedral spaces have also been calculated and found to be con- 
sistent with the WMAPl large-angle correlations (Aurich et al. 
I2005 i although for rather high values of the density parameter 
Qtot, in conflict with the current observations. 

The PDS model predicts pairs of matched antipodal cir- 
cles as a definite signature (Cornish et al. |1996 ). Roukema 
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( |2000al |2000bl ) was the first to apply the matched circle prin- 
ciple empirically to COBE data, whose poor angular resolution 
did not allow for a quantitative analysis. Next, different teams 
searched for such circles in the higher resolution WMAPl data, 
using various statistical indicators and massive computer calcu- 
lations. The team that first conceived the circle method obtained 
a negative result, which led them to reject the PDS model on a 
sub-horizon scale (Cornish et al. 2004; Key et al. 2007). A sec- 
ond team, using a different analysis, found six pairs of matched 
circles distributed in a dodecahedral pattern, with angular sizes 
consistent with the PDS model, but they provided no statisti- 
cal analysis so their claim remains unconvincing (Roukema et 
al. 2004). A third team (Aurich et al. 2005, 2006) performed a 
more careful search and analysis and found that the signal was 
considerably degraded by the integrated Sachs-Wolfe (ISW) and 
Doppler contributions to the CMB temperature fluctuations (see 
Eq. |22] below). They concluded that the PDS model could be 
neither confirmed nor rejected by the circle search on presently 
available data. Then ( 120061 ) also showed that even with fore- 
ground cleaned CMB maps, in which the most noisiest parts of 
the sky are masked, it is very hard to find matched circles. Thus, 
the topic remains controversial and the debate about the obser- 
vational pertinence of the PDS remains open. 

The recent release of the three-year WMAP results (Spergel 
et al. 2007; Page et al. 2007; Hinshaw et al. 12007 1 Jarosik et 
al. |2007 . hereafter WMAP3 papers) strengthens the evidence 
for weak large-angle correlations. The original WMAPl anal- 
ysis showed unusually weak CMB temperature correlations on 
angular scales greater than 60°, at a confidence level of 99.85% 
(Spergel et al. I2003I ). Still looking at scales greater than 60°, 
an independent analysis finds the correlations to be weak at the 
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99.91% level for the WMAPl dat a but th en at the 99.97% level 
for the WMAP3 data (Copi et al. '2006'), in both cases for the 
cut sky (which excludes local contamination from the Milky 
way). This strengthened evidence for weak large-angle correla- 
tions motivates continued interest in the PDS and related finite 
universe models. 

A multi-connected space is the quotient E - XfY of a sim- 
ply connected space X (the universal cover of E) under the ac- 
tion of a group r of symmetries of E (the holonomy group). For 
£■ to be a manifold (rather than an orbifold or something more 
complicated still), the holonomy group X must be discrete and 
fixed point free (see, e.g., Lachieze-Rey & Luminet l 19951 1. The 
universal cover X corresponds to a homogeneous spatial section 
of a standard Friedmann-Lemaitre spacetime. In the case of the 
PDS, the universal cover is the three-sphere S^, and the finite 
holonomy group T is generated by two elements, represented by 
two special orthogonal matrices. 

This paper evaluates further the predictions of a PDS cosmo- 
logical model and compares them with new CMB observations. 
This requires knowing the eigenmodes of the Laplacian on the 
PDS, i.e., the solution of its eigenvalue Helmoltz equation. Up 
until recently, only numerical solutions were available, and for 
the first modes only (corresponding to low wavenumbers k). For 
instance, the original calculations were purely numerical, and 
only up to A; = 24 (Luminet et al. '2003). Later work extended 
the calculation to ^ = 230 (Aurich et al. 2005). The present work 
is based on new analytical calculations of the eigenmodes of the 
PDS (Lachieze-Rey & Caillerie 2005), which are consistent with 
Bellon's recent work (2006). They allowed us to compute modes 
for very high values of k, much beyond k - 230. However, 
storage capacity limited the results presented in this article to 
k - 230 as the number of modes goes as k^^.^^, because a mode 
k is defined by the (k+ 1)^ coeflicients of its projection onto the 
modes of . If the aim is to calculate the power spectrum only 
(without simulating temperature maps), we can go further, up to 
k - 3) 000, using a theorem of Gundermann ( 2005J . which was 
also conjectured independently by Aurich et al. (|20051 l. 

The eigenmodes of the multi-connected space E lift to T- 
invariant modes of the universal cover X. Conversely, each T- 
invariant mode of X projects down to a mode of E. Thus, we 
may safely visualize - and compute - the modes of E as the T- 
invariant modes of X. They are the solutions ^ of the Helmholtz 
equation: 



(1) 



where A^ is the Laplacian on X and Ak is the eigenvalue asso- 
ciated to the integer wavenumber k. The invariance condition 
under F reads 



(2) 



for all symmetries g e F and all points x eX.\i follows that each 
eigenvalue for a multi-connected space E is also an eigenvalue 
for the universal cover X. When X - S^, the eigenvalues take 
the form Ak - -k(k + 2), indexed by integer wavenumbers k eN, 
with multiplicity (k + 1)^. The eigenvalues of the PDS form a 
subset of this set, given explicitly in Ikeda (1995). 

Given a model for gravitational instabilities, the statistical 
distribution of the CMB temperature fluctuations in the PDS 
model depends on the PDS's eigenmodes. Applying the method 
of Lachieze-Rey & Caillerie (2005), we construct the modes up 
to wavenumber A:n,ax - 230. From them, we calculate the implied 
realizations of the CMB temperature. This lets us reach a reso- 
lution in the temperature fluctuations corresponding to the (cur- 



vature dependent) angular wavenumber /"max ~ ^max V^^tot - 1 ~ 
30. 

In Sect. |3] we show typical maps corresponding to realiza- 
tions with angular resolution 6 ~ 6°. In these maps, we exhibit 
the expected presence of matched "circles in the sky" (Cornish 
etaL [T998] l. 

From the realizations of the temperature distributions, we 
have generated a statistical set of realizations of the eigenmodes 
distribution. This allowed us to estimate expectation values and 
statistics for the characteristics of the CMB temperature distri- 
butions, in particular for the angular power spectrum coefficients 
C(, up to the limit £,mx- 

For the modes with lowest values of { (up to 24), we confirm 
the previous calculations of Luminet et al. (120031 ). We extend 
them over a wider range of the spectrum and compare them with 
the WMAP3 data. 



2. Eigenmodes of the Laplacian 

We recall briefly the main results obtained in Lachieze-Rey & 
Caillerie (2005 ), summarizing the key elements and including 
practical details concerning the eigenmodes of the PDS, i.e., the 
solutions of Eqs. ([TJ and (|2]) for the PDS. 

The solutions of the Helmholz equation ([TJ for the 3-sphere 
S ^ form a vector space 'V, which is the direct sum of subspaces 



(3) 



where denotes the space of eigenmodes of with eigen- 
value Alt - -k(k + 2). Similarly, the vector space of eigenmodes 
of the PDS splits as the direct sum 



PDS 



PDS 



(4) 



k=a 



Each is a (possibly empty) subspace of the corresponding 

'Vk- The dimension of 'yPDS jg jj^g multiplicity of the eigenvalue 
Ak for the PDS, which is always an integer multiple of ^ -H 1 
and is non-zero only if k is even. Ikeda d 19951 1 has calculated the 
multiplicities explicitly. 

Thanks to the splitting (HJi, we will compute the eigenmodes 
on each 'yPDS separately. For convenience, we will express the 
modes of 'yPDS j-giative to the basis of 'Vk- 

First, we consider the parabolic basis of 'Vk, which was first 
introduced in Bander & Itzykson ( 119661 ) by group theoretical ar- 
guments. It corresponds to the following set of functions Tk-mi.m2, 
defined for A;, mi, m2 e N and -k/2 < {m\,m2) < k/2 by: 



a{cos X e''^)\^mx e'*)"'P™''(cos Ix) 



(5) 



where / :— m\ + OT2, m :— 
coefficient 



OTi, d - kjl - m2 and the 



k+1 (kl2 + m2)\(kl2-m2)\ 
27r2 ^{kl2 + mi)\{kl2-nn)\ 

is computed from normalization requirements. As T s {Ta;™, 
is a basis of 'V^, it generates 'V^, so each solution of dTJ for the 
PDS can be decomposed on this set as: 



(6) 
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where the fk;m,,m2 complex numbers. 

In order to be an eigenmode of the PDS, if)^^^ has to be 
invariant under F, the holonomy group of the PDS. To be so, it 
is sufficient that it be invariant under two generators, which we 
denote andg2: 

il>^^^(x)^^^^(g,(x)) and ^PDS(x) = ^™S(^^(_^)) 
for all X e 5 ^, or equivalently 

and (8) 

where the operator Rg expresses the natural left action of the 
SO(4) rotation g defined as {Rgij/){x) = t//{g^^x). 

To solve this system of equations, taking advantage of their 
nice properties under rotations on the sphere, we introduced the 
Wigner D-functions. This set of functions naturally acts on the 
Lie group SO(3) via the following action; 



D4,,„, : SO(3) 



< mi \Rg\j, 1112 > 



(9) 



The vectors | jm >, - j < m < j, form an orthonormal basis for 
'H^, the 2j + 1 dimensional irreducible representation of SU(2). 
When j is integer, this is also an irreducible unitary representa- 
tion of S0(3), and the | jm > can be taken as the usual spherical 
harmonics Yj,„. This implies 



jm'- 



(10) 



In Lachieze-Rey & Caillerie (2005'), we also derived a useful 
identity describing the action of these functions on the product 
of two elements of SO(3), 



(11) 



which is just the resolution of the identity in their definition 

The isometry between and SU(2), as manifolds, allows 
us to link these functions with those of the parabolic basis. Any 
point of is identified to an element of SU(2) by the following 
relation: 



x = {x,e,4>) 



Mf = 



SU(2) 

i smx e^'^ cosx ^ 



(12) 



Note that one finds other phase choices for this identification in 
the literature. 

On the other hand, there is a group isomorphism between 
SU(2)/Z2 and SO(3), where Z2 refers to the multiplicative group 
{+id, -id}. Thus, each element u of SU(2) defines a rotation gu 
in SO(3). In practice, a rotation of SO(3) is parameterized by 
its Euler angles a, /3, y. Taking into account the identification 
above, the correspondence takes the form: 



SU(2) 
with 

.4 . 



SO(3) 



a - y 



(14) 



This lets us lift the Wigner D-functions from SO(3) to SU(2) 
or, equivalently, to S^. Their explicit expression (see for instance 



Edmonds 1 19601 1 shows the identification, up to a constant, with 
the previous functions Tt,mum2'- 



27r2 
k+1 



Tk\m] ,1112 (^)j 



(15) 



with -kl2 < m\,m2 < k/2. 

We can now write Eq. (|7J explicitly, by expanding modes on 
the parabolic basis. This gives, for any x in S^, 

^ //:;mi ,m2 -^A:;mi ,m2 (■^) ~ ^ ] fk\m\,m2 T^k\m\,m2^S ■^) (1^) 



where g is one of the two generators {gi or ^2) of the PDS (or, 
if desired, of some other spherical space). Using property ( fTTT i. 
we obtain a simple relation between the coefficients of the eigen- 
mode s: 



fk;mi,m2 — ^^fk;mi,m I^nim2^S )• 



(17) 



Note that this equation is independent of the index nii, so we can 
further split each vector subspace as a direct sum 



(18) 



More practically, this corresponds to writing the eigenmode 

fk;n 



as 



fk;n 



fkm2 ^m(,m\ 7 



(19) 



where niQ is some arbitrary integer between -k/2 and k/2. This 
implies that the multiplicity of an eigenmode of a spherical space 
is proportional to k + 1. Thus, we may rewrite Eq. ( [TT] ) as an 
eigenvalue equation 



Dk.gFk - Fk 



(20) 



where Fk is the vector with coordinates fkm2 and Dk^g is the ma- 
trix whose general term is Z)f/m2(^ ')■ This gives a system of two 
eigenvalue equations to be solved simultaneously. In Lachieze- 
Rey & Caillerie ( 2005 ) , we have shown that when we choose a 
basis of the three-sphere where one generator, say gi, is diagonal 
(which is always possible), this imposes the condition 



m2 = (mod «) 



(21) 



where the integer n equals 5 in the case of the PDS. In this case, 
finding the eigenmodes of the Laplacian corresponds to finding 
a solution of ( l20t for g = gi and where we consider only indices 
1112 satisfying the condition i2T[ . 



3. Temperature maps 

We now restrict our study to cosmological models whose spatial 
section is a PDS, and to values of iltot in the range favored by the 
WMAPl and WMAP3 observations, namely Qtot = 1 .02 + 0.02. 



(13) 3.1. Maps 



Matter fluctuations at recombination are assumed to follow a 
Gaussian random distribution with a Harrison-Zel'dovich power 
spectrum with P(k) cc k^. Thus, a realization is given by a 
Gaussian random distribution of the modes calculated in the 
previous section, constrained to have the desired power spec- 
trum law. However, as first emphasized by Roukema (2000;^, 
the power spectrum approaching the injectivity diameter and the 
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out-diameter of the fundamental domain may well be different 
from the theoretical expectation, since these scales represent 
the physical size of the whole universe, and the observational 
arguments for a A:' spectrum at these scales are only valid by 
assuming simple connectedness. This can be considered as a 
caveat for the interpretation of the following results. 

Such a distribution of matter fluctuations generates a temper- 
ature distribution on the CMB that results from different physical 
effects. If we subtract foreground contamination, it will mainly 
be generated by the ordinary Sachs-Wolfe (OSW) effect at large 
scales, resulting from the the energy exchanges between the 
CMB photons and the time-varying gravitational fields on the 
last scattering surface (LSS). At smaller scales, Doppler oscilla- 
tions, which arise from the acoustic motion of the baryon-photon 
fluid, are also important, as well as the OSW eff'ect. The ISW ef- 
fect, important at larger scales, has the same physical origin as 
the OSW effect but is integrated along the line of sight rather 
than on the LSS. This is summarized in the Sachs-Wolfe for- 
mula, which gives the temperature fluctuations in a given direc- 
tion h as 

6T Il6p \ n" . . 

— («)= T— + 3> ('7LSs)-«.v.(/7LSs)+ (a) + >F)d77 (22) 
T \4 p I J^^^^ 

where the quantities <i> and *F are the usual Bardeen potentials, 
and Vg is the velocity within the electron fluid; overdots denote 
time derivatives. The first terms represent the Sachs-Wolfe and 
Doppler contributions, evaluated at the LSS. The last term is 
the ISW effect. This formula is independent of the spatial topol- 
ogy, and is valid in the limit of an infinitely thin LSS, neglecting 
reionization. 

The temperature distribution is calculated with a CMBFast- 
like software developed by one of u£], under the form of temper- 
ature fluctuation maps at the LSS. One such realization is shown 
in Fig. [H where the modes up to - 230 give an angular res- 
olution of about 6° (i.e. roughly comparable to the resolution 
of COBE map), thus without as fine details as in WMAP data. 
However, this suffices for a study of topological effects, which 
are dominant at larger scales. 

Such maps are the starting point for topological analysis: 
firstly, for noise analysis in the search for matched circle pairs, 
as described in Sect. 13.21 secondly, through their decompositions 
into spherical harmonics, which predict the power spectrum, as 
described in Sect. |4] In these two ways, the maps allow direct 
comparison between observational data and theory. 

3.2. Circles in the sky 

A multi-connected space can be seen as a cell (called the fun- 
damental domain), copies of which tile the universal cover If 
the radius of the LSS is greater than the typical radius of the 
cell, the LSS wraps all the way around the universe and inter- 
sects itself along circles. Each circle of self-intersection appears 
to the observer as two different circles on different parts of the 
sky, but with the same OSW components in their temperature 
ffuctuations, because the two different circles on the sky are re- 
ally the same circle in space. If the LSS is not too much bigger 
than the fundamental cell, each circle pair lies in the planes of 
two matching faces of the fundamental cell. Figure |2] shows the 
intersection of the various translates of the LSS in the universal 
cover, as seen by an observer sitting inside one of them. 



' A. Riazuelo developed the program CMBSlow to take into account 
numerous fine efi'ects, in particular topological ones. 




Fig. 1. Temperature map for a Poincare dodecahedral space with 
Otot - L02, - 0.27 and h - 0.70 (using modes up to 

k - 230 for a resolution of 6°). 




Fig. 2. The last scattering surface seen from outside in the uni- 
versal covering space of the Poincare dodecahedral space with 
Qtot - L02, Qmat = 0.27 and h - 0.70 (using modes up to 
k - 230 for a resolution of 6°). Since the volume of the physical 
space is about 80% of the volume of the last scattering surface, 
the latter intersects itself along six pairs of matching circles. 



These circles are generated by a pure Sachs-Wolfe effect; in 
reality additional contributions to the CMB temperature fluctua- 
tions (Doppler and ISW effects) blur the topological signal. Two 
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teams have carefully analyzed the blurring in the case of the 
Poincare dodecahedral space: one team finds it strong enough 
to hide matching circles (Aurich et al. 2005, 2006), while the 
other team reaches the opposite conclusion (Cornish et al. 120041 
Shapiro Key et al. i2007 ). By contrast, in the case of a small 3- 
torus universe, everyone agrees that blurring could not hide cir- 
cles. The blurring's effectiveness varies by topology, because the 
relative strengths of the various effects (OSW, ISW, Doppler) de- 
pend on the modes of the underlying space. 

4. Power spectrum 

Like any function defined on the sphere, the CMB temperature 
fluctuations can be decomposed into spherical harmonics 

ST ' 

T^"^ = Z Z ^'"'^"^ ^2^^ 

/=0 m=-l 

where the unit vector h e represents a point on the sky. 

For a statistically homogeneous and isotropic distribution of 
matter in the universe, the power spectrum 

1 ^ 

m=-l 

contains all relevant information about the temperature fluctua- 
tions. Such a power spectrum has been calculated for the CMB 
temperature fluctuations as measured by WMAP. We compare it 
here with the predictions of our model. 

From our realizations of the predicted temperature maps, 
we have calculated power spectra. Assuming a present day 
matter density Qmat - 0.27 and a reduced Hubble parameter 
h - 0.70 (Spergel et al., 120071 ), we have varied Qtot within the 
range 1.015-1.025, where we get low-/ power compatible with 
WMAP3 data. For illustrative purposes. Fig. [3] shows our esti- 
mated spectra for the PDS, using three values of fltot within the 
above range, compared with the WMAP3 data taken from Fig. 2 
of Spergel et al. ( 120071) . We have normalized the curves by insist- 
ing that they approach the concordance model curve for high {, 
because topology is significant only for low f. We can see that 
the PDS model agrees well with observations for all three val- 
ues of Qtot, up to the cutoff /"max ~ ^max V^tot - 1 • Sincc our 
numerical computations go only to wavenumber k^^^ ~ 230, the 
resulting power spectrum is reliable only up to ^ ~ 25 - 30, de- 
pending on Qtot. The oscillatory pattern in the C/ is a prediction 
of the Poincare dodecahedral space model due to multiplicities 
by which the vibrational modes are weighted in the mode sum. 
The effect is significant only for the low vibrational modes; it 
is washed out at sufficiently large I (/ > 30) so that the asymp- 
totic behaviour of the power spectrum recovers the curve for the 
simply connected . 

Figure |4] shows the optimal fit at Qtot = 1.018 for which 
the quadrupole suppression is maximal, a value perfectly con- 
sistent with the WMAP team reported Qtot = 1.02 ± 0.02. In 
order to compare the predictions of simply connected and multi- 
connected topologies, we also compare our best fit PDS spec- 
trum with the standard best fit "concordance model" spectrum. 

5. Conclusion 

Recent analytical eigenmode calculations for the Poincare do- 
decahedral space allowed us to simulate CMB temperature fluc- 
tuation maps more accurately. We confirmed the correctness of 
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Fig. 3. Comparative power spectra (in fjK^) as a function of the 
multipole f for WMAP3 (errorbars) and PDS for three values 
of Qtot (Qtot = 1.015 for the dashed curve, Qtot - 1.02 for the 
solid curve and Qtot = 1.025 for the dotted curve), assuming 
Qnjat = 0.27 and h - 0.70. Here we calculate the modes up to 
k = 3000 using the conjecture of Aurich et al. (i2005j proved by 
Gunderman ( |^05] l. 
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Fig. 4. Comparative power spectra (in yuK^) as a function of the 
multipole { for WMAP3 (errorbars), "concordance model" (dot- 
dashed curve) and PDS (solid curve) for Qtot = 1.018, Qmat = 
0.27 and h = 0.70. 



the maps by verifying the presence of the expected circles-in- 
the-sky in the OSW-only maps. 

Using a random set of Gaussian realizations of the matter 
fluctuations, we have calculated the predicted power spectrum 
of the CMB temperature fluctuations and the two-point temper- 
ature correlation function. Our results for the lowest modes con- 
firm the numerical estimates of Luminet et al. ( 2003 ), while use 
of higher PDS modes let us estimate the multipoles up to ^ ~ 30. 
We have obtained an excellent fit with the WMAP data, imply- 
ing that the PDS cosmological model remains a good candidate 
for explaining the angular spectrum, even though the negative re- 
sults of matching circle searches remain a topic of debate (Key 
et al. 2007, Aurich et al. 2006, Then 2006). 

Clearly the power spectrum alone cannot confirm a multi- 
connected cosmological model. Although the PDS model fits 
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the WMAP3 power spectrum better than the standard flat in- 
finite model does, alternative explanations may still be found, 
the simplest one being an intrinsically non-scale invariant spec- 
trum. Thus, one must find complementary checks of the topolog- 
ical hypothesis. The off-diagonal terms of the correlation matrix 
provide one possibility. Unfortunately, they are not easy to de- 
rive and are difficult to interpret statistically. On the other hand, 
polarization eff'ects provide an additional tool for investigating 
space topology (Riazuelo et al. 120061 1. 
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